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D . 
pL( . Abstract 

We consider the SUSY-model neutralino as a dominant Dark Matter 
pq ■ particle in the galactic halo and investigate some general issues of direct 

^ \ DM searches via elastic neutralino-nucleus scattering. On the basis of con- 

ventional assumptions about the nuclear and nucleon structure, without 
p^ ■ referring to a specific SUSY-model, we prove that it is impossible in prin- 

ciple to extract more than three constraints on fundamental SUSY-model 
parameters from the direct Dark matter searches. Three types of Dark 



o 



Q^ ■ matter detectors probing different groups of parameters are recognized. 
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Many physicists at present believe that more than 90% of the total universe 
mass should be in the form of non-luminous non-baryonic dark matter (DM). 
Most of our galaxy mass should be in a dark halo as well. According to the 
commonly accepted halo model the DM particles have a Maxwellian velocity 
distribution in the galactic frame, a mass density in the Solar system of about 
p ~ 0.3 GeV-cm~^ and come to the Earth with mean velocities v ~ 320 km/sec, 
producing a substantial flux $ = p ■ v/M ($ > 10^ cm~^ sec~^ for a particle 
mass M ~ 1 GeV). Thus the idea of detecting DM particles directly, for instance 
through the elastic scattering off nuclei in detector target material looks realistic. 

Dark Matter in the universe consists of two main components, which play 
different cosmological roles. These are "hot" (HDM) and "cold" (CDM) DM 
components which decoupled from the thermal equilibrium in the early universe 
at relativistic and non-relativistic temperatures, respectively. The analysis of the 
COBE satellite data |T[] on anisotropy in the cosmic background radiation and 
the theory of the formation of large scale structure in the universe lead to the 
following DM composition DM = 0.7(CDM) + 0.3(HDM) g], ||. It is commonly 
believed that the CDM is represented by weakly interacting massive particles 
(WIMP) and the axion, while the HDM consists of light neutrinos. We do not 
discuss the axion as a DM candidate, concentrating on the pure WIMP variant 
of the CDM. 

The neutralino (x) is a SUSY-model version of WIMP. Presently it is the most 
favorable CDM candidate, x is a Majorana (x^ = x) particle with spin | and 
usually assumed to be the lightest SUSY-particle (LSP). 

The problem of direct detection of the DM neutralino x via elastic scattering 
off nuclei has been considered by many authors and remains a field of great 



experimental and theoretical activity p|-p!3|. 

One can meet in the literature the statement that the neutralino prefers to 
scatter off spin-non-zero (Jt^O) nuclei due to its Majorana nature. However, we 



have shown ||13[ that for sufficiently heavy nuclei nuclear spin is a non-dominant 
factor in the neutralino-nucleus scattering. Thus, one has equal chances to detect 
the DM signal using either spin- zero or spin-non-zero target nuclei with atomic 
weight A > 50. 

However, Jt^O and J=0 - nuclei provide different information about SUSY- 
model parameters. Therefore experiments with Jt^O and with J=0 - nuclei are 
both important from the point of view of SUSY-phenomenology. 

In the present paper we envisage the question of how many independent con- 
straints on fundamental SUSY-model parameters can be extracted from direct 
DM neutralino searches. On the basis of conventional assumptions we shall prove 
that there are just three independent constraints for special groups of parameters. 



In this respect we can distinguish three types of DM detectors with different 
sensitivity to the groups of fundamental parameters. These are detectors with 
non-zero-spin nuclei of neutron-like ('^He, ^^Si, ^^Ge,...) and proton-like (^^F, 
35q^ 205'jj^ -j shell-model structure, and spin-zero nuclei (^^Ge,...). 

The DM event is elastic neutralino-nucleus scattering causing a nuclear recoil 
detected by detector. The event rate per unit mass of target material depends on 
the density and velocity of DM in the solar vicinity and the cross section o"e/(x^) of 
the neutralino-nucleus (x^) elastic scattering. One can calculate aeiix^) starting 
from the neutralino-quark effective Lagrangian. It can be written in the general 
form as 

^eff = II A ■ Xl^^l5X ■ qi^im + Yr~ -^1-XX- QQ, (1) 



M, 



w 

where terms with the vector and pseudoscalar quark currents are omitted being 
negligible in the case of non-relativistic DM neutralinos with typical velocities 
f^ ~ 10~^c. 

The coefficients Aq, Sq are effective phenomenological parameters completely 
parametrizing low-energy neutralino-quark interactions. They can be calculated 
in a specific SUSY-model. Furtheron we will consider the most general properties 
of the x-^ scattering starting from the effective Lagrangian (|I|) without referring 
to specific SUSY-models. To calculate the elastic neutralino-nucleus cross sec- 
tion aei{x^) we follow the standard procedure of sequentially averaging the x-Q 
interactions over the nucleon and then the nuclear structure. 

The first and the second terms in L^ff (^ averaged over the nucleon states give 
spin-dependent and spin- independent matrix elements Mgd and Msi, respectively. 
For the spin-dependent matrix element we have 0,0: 

Mfi")=44s;(„)->Vp(.), (2) 

where S^ and S'p(n) are the neutralino and proton (neutron) spin-operators. The 
factor Wp(„) depending on the parameters of the effective Lagrangian (0) has the 
form 

>Vp(n) = E A^Q^ (3) 

qep{n) 

Aq are the fractions of the nucleon spin carried by the quark q. The definition is 

< p{n)\qYl5q\p{n) >= 2S;(„)Ag, (4) 

where S'p/^N = (0, 5^(n)) is the 4-spin of the nucleon. The parameters Aq for the 
proton can be extracted from the EMC |[T^ and hyperon data [|T5|: 



Am = 0.77 ±0.08, Arf = -0.49 ± 0.08, As = -0.15 ± 0.08. (5) 



The relevant values for the neutron can be found from (|^) by isospin symmetry 
substitution u ^ d,d ^ u. 

The spin- independent matrix element M^j has the form [0!@]Q- 

M,, = ^-xX-^^-V. (6) 

Again as in the spin-dependent matrix element (^ we separate the factor 

v = Y.f,-s„ (7) 

where the parameters fg are defined as follows: 

< pin)\mg ■ qq\pin) >= fg ■ Mp(„)^^. (8) 

The values extracted from the data are |jT6|,[|l^: 

/„ = 0.022, /, = 0.034, fs = 0.38, 
2 

Jc^Jb^Jt^ T^ \^ Ju Jd Js)- 

The next step is the calculation of nuclear matrix elements of scalar and axial- 
vector nucleon operators ^ ■ \l/ and ^7^75\l/, respectively for the nuclear state 
\A >. The axial- vector operator corresponding in a non-relativistic limit to the 
nucleon spin-operator appears in eqn. (0). We can parametrize these nuclear 
matrix elements in the most general form as follows 



<A|Mp(„)M|A> = MAFs^{q^)AA, (9) 



<A|4,(„)|A> = A,(„)Ffi")(g2)<A|J>>. 



Here J is the nuclear spin. The form factors Fsfi,si{q^) take into account effects 
of finite momentum transfer q and obey the normalization conditions: 



sd,si 



(0) = 1. (10) 



For nuclei with atomic weight A > 50 nuclear structure is essential since qmax > 
1/Ra [Ra is the nuclear radius) [1^. In this case the form factors strongly deviate 
from the normalization value. 



^ When this paper had been completed we received a paper ref . [|ll| with more refined treat- 
ment of the spin- independent matrix element. 



The phenomenological parameter Xp(n) measures the contribution of the pro- 



ton (neutron) group to the nuclear spin (See [18| for these nomenclature.). Ac- 
cording to the conventional shell model picture, the odd-group contribution dom- 
inates for even-odd nuclei. For odd-proton group nuclei (^Li, ^^F, ^'^^Tl,...) 
Xp >> \n and for odd-neutron group nuclei (^He, ^^Si, ^^Ge,...) A„ >> Xp. 

In the case of odd-odd nuclei (^°B, ^^N,...) proton and neutron groups give 
comparable contributions to a nuclear spin and A„ ~ Xp. 

The A-parameters should be calculated in some nuclear model. In the semi- 



empirical odd-group shell model [|18l they can be related to the nuclear magnetic 



moment, /i, as follows 








XJ = 


_ fi- 9^J 
9' - 9^' 



:ii: 



where g^ = 1(0) and g^ = 5.586(— 3.826) are orbital and spin proton (neutron) 
(yf-factors. 

Then one can extract A- values from the experimental data on nuclear magnetic 
moments^. 

For our purpose the numerical values of these parameters are irrelevant but 
for completeness we give them in the Appendix for nuclei of interest in the DM 
search. 

Making standard transformations of the above-written matrix elements and 
averaging over the Maxwellian distribution of neutralino velocities, we finally 
arrive at the formula for the event rate of the elastic neutralino-nucleus scattering 
in the detector per day per unit mass of the target material: 



R 



K + an-K + ao-V' 



events 
kg -day' 



:i2) 



The parameters Oj depend on properties of the target nucleus as well as on the 
mass density and the average velocity of DM particles in the solar vicinity. In 
the Appendix we give the values of these parameters for different nuclei within 
specific model assumptions. The above-defined quantities W and V contain all 
dependence on the parameters of the effective Lagrangian (|l]) and do not depend 
on nuclear properties. This factorization is a central point of the present paper. 
It follows from eqn. (|12]) that measuring the event rate R we can study just 
three special combinations of fundamental parameters Wp^n and V defined in 
eqn. (^), (|^). This is the only information about fundamental parameters ac- 
cessible in DM search experiments. R is a linear combination of the quantities 

^Another approach based on the theory of finite Fermi systems is described in ||l^ . 



W^p,„ and V^. To extract experimental limitations for each of them one should 
search for DM with different target nuclei. We can distinguish three types of DM 
detectors by their sensitivity to Wp^n and V. These are detectors built of: 

1. spin-non-zero target nuclei with an 'odd-proton group' probing a linear 
combination 

ap-W'p+ao-V'<R% 

and giving an experimental limit R^Jp] 

2. spin-non-zero target-nuclei with an 'odd-neutron group' probing another 
linear combination 

an-Wl + ao-V'<RfJp 

and giving an experimental limit R^Jp] 

3. spin-zero target nuclei sensitive only to the scalar part of the neutralino- 
nucleus interaction 

K ^ i ''exp 

giving an experimental limit R^J^X- 

In conclusion we would like to stress the following. To extract all information 
about SUSY-model parameters possible from direct DM search one should have 
three established types of DM detectors. No other information can be obtained 
from the direct DM search experiments. Different detectors can only improve the 
data on the three above-defined groups of SUSY-model parameters. 

A detailed discussion on the extractable information about SUSY model pa- 
rameters will be given elsewhere. 



1 Appendix 

Below we give complete formulae for calculation of the nuclear structure param- 
eters entering in the definition of the event rate R in eqn. (|1^) within commonly 
accepted model assumptions. 

Assuming a gaussian parametrization |10| for the form factors Fsd,si{<f) in 
eqns. (^ and averaging over the Maxwellian distribution of DM particle velocities 
we get 

ap(„) = 5.8-10i°-A^(„)J(J+l)-C(r.pJ-P, (13) 

a, = 1.44-10i°-(^)\(r,,,,,,)-P, (14) 

iviw 



where the DM factor is 



V 



^M^Ma 



P 



(M^ + M^)2J L.3Ge1/ ■ cm^sJ V^2Qkm/s 



< \ve\ > 



(15) 



Here M^,Ma are the DM neutrahno and target nucleus masses. The coherence 
loss factor EO] is defined by 



^,^,^aW3/ _exp(-^erf(^ ,16) 



b V VT+b erf(l) 

where 

= —a r 



9 (M^ + M^)2- 

Here o"^ is the dispersion of the Maxwellian neutralino velocity distribution a = 
0.9 ■ 10~^. r^pin and r charge are the rms spin and charge radius of the nucleus A. 



To estimate r charge one can use the parametrization ||2T| : 

reharge = (0.3 + 0.89M'f) fm. (17) 

The values of the ratio Tspm/rcharge can be estimated in the harmonic oscillator 
shell model 0. 

Numerical values of the above defined parameters for nuclei of interest in 
direct DM searches are given in the Table. 
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Table: Parameters for event rate calculations. 



Isotope 


J 


H 


K 


''^spin 


Isotope 


J 


H 


K 


"^ spin 


''"charge 


''"charge 


iH 


1/2 


1.0 





1.00 


99Ru 


5/2 





0.0045 


1.19 


^Hc 


1/2 





1.2373 


1.00 


loiRu 


5/2 





0.0056 


1.19 


^Li 


3/2 


0.1096 





1.17 


107 Ag 


1/2 


0.0720 





1.06 


^Be 


3/2 





0.0768 


1.12 


109 Ag 


1/2 


0.0760 





1.06 


i^B 


3/2 


0.0299 





1.09 


ii^Cd 


1/2 





0.0960 


1.17 


15N 


1/2 


0.1160 





1.06 


"^Cd 


1/2 





0.1053 


1.17 


1-^0 


5/2 





0.0391 


1.25 


ii^Sn 


1/2 





0.2307 


1.16 


19p 


1/2 


0.8627 





1.21 


"7Sn 


1/2 





0.2733 


1.16 


23Na 


3/2 


0.0109 





1.16 


i^^Sb 


5/2 


0.0057 





1.15 


^^Al 


5/2 


0.0099 





1.13 


i23Sb 


7/2 


0.0035 





1.15 


^^Si 


1/2 





0.0840 


1.12 


127J 


5/2 


0.0026 





1.15 


31p 


1/2 


0.0760 





1.11 


i29Xe 


1/2 





0.1653 


1.14 


35C1 


3/2 


0.0096 





1.10 


i^iXe 


3/2 





0.0147 


1.14 


47Ti 


5/2 





0.0067 


1.20 


i^^Cs 


7/2 


0.0033 





1.14 


49Ti 


7/2 





0.0068 


1.20 


i39La 


7/2 


0.0020 





1.13 


51V 


7/2 


0.0106 





1.19 


i^^Gd 


3/2 





0.0021 


1.22 


55Mn 


5/2 


0.0069 





1.17 


i^^Gd 


3/2 





0.0035 


1.22 


59 Co 


7/2 


0.0049 





1.15 


183^ 


1/2 





0.0040 


1.19 


e^Zn 


5/2 





0.0083 


1.13 


i9ilr 


3/2 


0.0387 





1.08 


69Ga 


3/2 


0.0056 





1.13 


193 J, 


3/2 


0.0379 





1.08 


71 Ga 


3/2 


0.0237 





1.13 


199Hg 


1/2 





0.0693 


1.17 


7^Ge 


9/2 





0.0026 


1.25 


20lHg 


3/2 





0.0096 


1.17 


793, 


3/2 


0.0077 





1.11 


203r|nj 


1/2 


0.2400 





1.07 


«iBr 


3/2 


0.0125 





1.11 


205'jj 


1/2 


0.2467 





1.07 


91 2r 


5/2 





0.0186 


1.21 


2U7pb 


1/2 





0.0960 


1.17 


93Nb 


9/2 


0.0065 





1.20 


209Bi 


9/2 


0.0002 





1.16 



